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THE  INADEQUACIES  OF  DAMAGE  ENERGY  AS  A  MEASURE 
OF  DISPLACEMENT  DAMAGE 


INTRODUCTION 

One  of  the  few  simple  measures  of  radiation  damage  that  relates  damage  in  crystals  at  the 
microscopic  level  to  the  gross  behavior  of  the  material  is  the  number  of  Frenkel  pairs  produced 
by  a  given  irradiation.  Even  though  a  significant  fraction  of  the  pairs  initially  produced  quickly 
recombine,  the  expectation  is  that  the  number  of  vacancies  produced  by  an  irradiation  is  a  rough 
indication  of  the  extent  of  permanent  damage  to  the  material.  Currently  there  is  interest  in  the 
calculation  of  displacement  damage  because  of  the  apparent  relationship  between  the  number  of 
displacements  produced  by  an  impinging  ion  and  the  damage  suffered  by  semiconductor  materials 
and  devices.’*^ 

The  underpinning  of  calculations  of  radiation  damage  is  a  model  of  the  isolated 
displacement  event.  Suppose  an  ion  of  charge  Z^e  and  mass  is  moving  with  energy  E  in  a 
material  and  collides  with  a  lattice  atom  of  charge  Z^e  and  mass  A^,  transferring  an  energy  T  to 
the  bound  atom.  The  conditions  under  which  a  struck  atom  leaves  its  native  site  depends  on  the 
momentum  of  the  impinging  ion,  on  the  crystal  structure  of  the  material  and  on  the  form  of  the 
potential  that  holds  the  material  together.  A  number  of  general  discussions  of  displacement 
effects  are  available  (Dienes  and  Vineyard,^  Sosin  and  Bauer,^  Robinson,®  Lehmann,®  Gittis,^ 
Averback  and  Kirk^) 

To  delineate  our  displacement  model,  we  define  a  series  of  energy  parameters.  Let  Eg  be 
the  displacement  energy;  if  T  is  less  than  Eg,  then  the  struck  atom  stays  in  place  and  shares  its 
newly  gained  energy  with  the  lattice,  but  if  T>Eo,  then  the  struck  atom  is  displaced  and  precedes 
to  move  through  the  crystal.  (In  reality  the  displacement  energy  in  a  crystal  depends  on  direction, 
but  we  assume  an  isotropic  value.)  Further,  if  the  incoming  ion  transfers  an  energy  T>Eo,  but 
retains  an  energy  E-T<Eo,  then  it  replaces  the  struck  atom  at  the  lattice  site--no  vacancy  is 
produced.  If  more  than  one  species  is  present  in  the  calculation,  the  role  of  replacements  is  more 
complicated.) 

We  next  allow  a  binding  energy  Eg.  If  a  vacancy  is  produced,  the  displaced  atom  retains 
energy  T-Eg;  the  binding  energy  Eg  is  accounted  for  by  the  relaxation  of  the  lattice  atoms 
surrounding  the  vacancy.  We  assume  that  Eg<Eo  and  that  the  effect  of  the  binding  energy  is 
assessed  after  the  collision.  For  example,  if  T>Eo,  then  there  is  a  displacement  even  though  it 
may  be  the  case  that  T-Eg<Eo.  Finally,  we  presume  there  are  no  binding  losses  in  replacement 
collisions. 

Mmucript  approved  January  29,  1989. 
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We  will  follow  the  course  of  Robinson  and  Torrens’  in  their  binary-collision  simulation 
code  MARLOWE  by  defining  a  cutoff  energy  Ec;  ions  moving  with  energies  less  than  E^  are 
considered  to  have  stopped.  This  concept  is  useful  in  simulation  calculations  because  it  saves  the 
time  necessary  to  follow  the  numerous,  but  ineffectual,  slowly  moving  ions.  While  such  a  cutoff 
has  no  particular  benefit  for  analytical  displacement  calculations,  we  retain  it  so  as  to  ease 
comparisons  of  our  work  with  the  simulation  models.  We  will  use  the  values  Ec>  0,  Eg,  lEg. 

Some  of  an  ion’s  initial  energy  E,  a  portion  labeled  1(E),  will  go  toward  exciting  the 
electrons  in  the  solid.  Robinson  and  Torrens  categorize  the  bins  into  which  the  remaining  energy 
falls  as  follows:  binding  energy  loss  B(£),  as  discussed  above;  remaining  kinetic  energy,  the 
energy  retained  by  ions  moving  with  energies  less  than  E^;  and  subthreshold  loss,  energies  less 
than  Eg  acquired  by  (undisplaced)  lattice  atoms.  We  will  combine  the  last  two  quantities  into  one, 
which  we  call  the  thermal  loss  H(E).  Thus  we  have 

E  -  1(E) +  D(E)  -  1(E)  +  B(E)  +  H(E)  (1 

The  quantity  0(E)  is  usually  called  the  damage  energy,  which  by  definition  is  all  the  energy  not 
lost  to  inelastic  processess. 

Robonson  and  Oen^’  point  out  that  some  of  the  subthreshold  energy  is  lost  to  inelastic 
processes,  so  that  our  1(E)  should  be  somewhat  larger  and  H(E)  somewhat  smaller.  This  effect  has 
no  consequence  for  our  principal  results,  which  depend  on  the  direct  calculation  of  the  number  of 
displacements  without  using  the  damage  energy  as  an  intermediate  step.  Their  modification 
would  make  some  difference  in  our  estimates  of  damage  energy,  but,  as  discussed  later,  the 
concept  of  inelastic  losses  at  very  low  energies  is  of  dubious  value. 

In  the  early  calculations  of  displacement  damage^ ^  the  inelastic  loss  was  not  specifically 
included.  Instead,  the  partition  was  often  made  that  all  of  an  ion’s  energy  above  some  value  E, 
was  lost  to  the  electrons,  while  all  of  its  energy  below  this  value  went  into  creating  displacements 
by  means  of  elastic  collisions;  that  is,  the  damage  energy  was  E  for  E<E|  and  E]  for  E>E|.  As  the 
calculations  became  more  sophisticated,  the  inelastic  losses  were  included  throughout  the 
calculation.  The  loss  1(E)  to  the  electrons  could  be  directly  calculated,  but  the  remaining  energy 
was  still  called  the  damage  energy.  It  is  not  obvious,  given  the  simultaneous  inelastic  losses  as  the 
displacements  are  produced,  that  there  is  a  direct  relation  between  the  damage  energy  and  the 
number  of  displacements.  The  examination  of  this  point  is  the  main  goal  of  this  report. 

Several  other  assumptions  should  be  pointed  out.  For  a  lattice  atom  that  is  given  energy 
T<Ec,  we  assume  that  all  of  the  energy  T  goes  into  the  thermal  energy  bin  H(E).  For  bound  or 
captured  ions  with  kinetic  energies  less  than  Eg,  we  assume  all  of  the  energy  is  lost  thermally . 
even  if  these  energies  are  greater  than  Eg. 

Turning  specifically  to  displacement  mechanisms,  as  opposed  to  energy  loss  mechan'sms, 
we  note  that  an  ion  moving  in  a  crystal  with  energy  less  than  Eg  can  produce  neither  replacements 
nor  vacancies.  If  the  ion  has  an  energy  between  Eg  and  2Eg,  it  cannot  produce  a  vacancy,  but 
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may  generate  a  replacement  collision.  Figure  1  summarizes  the  coUisional  possibilities  of  our 
model.  Starting  with  a  zero  cutoff  energy  in  fig.  la,  we  see  the  thermal  loss,  vacancy  production 
and  replacement  regions  of  the  E-T,T  plane.  For  an  ion  of  energy  E  (E>S.4Eo  in  this  example), 
all  of  its  possible  (elastic)  collisions  lie  on  the  dashed  4S  degree  line.  In  fig.  lb  we  impose  the 
cutoff  presume  that  all  ions  of  energy  less  than  the  cutoff  lose  their  energy  to  thermal 

processes;  specifically,  there  are  no  electronic  losses  in  the  cutoff  (C)  region.  We  see  that  the  only 
difference  in  Eq-O  and  Ec^Eq  calculations  is  in  the  relative  amounts  of  thermal  and  electronic 


V:  VACANCIES 
R;  REPLACEMENTS 
T:  THERMAL  LOSSES 
C:  CUTOFF  REGION 


Fig.  1  "  The  various  collision  possibilities  as  a  function  of  energy  transfered  and  energy 
retained;  shown  for  cutoff  energies  of  a)  E(^,  b)  Ep-Eo  and  c)  E.— 2Eo.  The  thermal  loss, 
replacement  producing,  vacancy  producing,  and  cutoff  regions  are  shown.  See  the  text  for  a 
detailed  discussion. 
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loss.  As  we  view  the  case  Ec*2Eo  in  fig.  Ic  we  note  a  further  shift  in  the  relative  thermal  and 
inelastic  losses.  A  cutoff  of  2Ed  also  reduces  the  estimated  number  of  replacements,  because  some 
ions  that  would  undergo  replacement  collisions  are  no  longer  followed. 

Almost  all  displacement  calculations  have  assumed  that  the  initial  ion,  of  energy  E>Eo,  is 
itself  a  displacee — produced,  for  example,  by  an  incident  neutron.  If  N(E)  is  the  average  number 
of  vacancies  resulting  from  an  ion  of  energy  E,  then  this  initial  displacement  assumption  amounts 
to  the  requirement  that 

f  0,  E<Eo 

N(E)  -  \  (2 

(.1,  Eq  <  E  <  2Eo  . 

We  find  it  more  natural  to  assume  the  cascade  starts  with  an  ion  of  energy  E  moving  freely  in  the 
crystal  and  to  calculate  the  number  of  secondary  vacancies  that  are  produced,  which  we  call 
Nv(E).  For  the  monatomic  substances  that  we  consider,  we  find— it  follows  from  the  integral 
equation— that  Nv(E)«0  for  E<2£o.  The  two  estimates  of  displacement  damage  are  related  by 
N(E)-Nv(E)+l  for  E>Eo. 

The  earlier  calculations  of  displacement  damage  involved  solving  an  integral  equation^  for 
N(E).  More  recently,  the  LNST  (Lindhard,  Nielsen,  Schraff  and  Thompsen)^^  results  for  the 
damage  energy  D(E)  (see  Sigmund^^  for  a  brief  history  of  their  approach)  have  been  widely  used 
to  estimate  displacement  damage  through  the  use  of  the  Khinchin- Pease  formula 

N(E)  -  «  EKE)  /  (2  Ed)  ,  (3 

where  /c  is  a  constant  of  the  order^^  of  0.8.  In  our  approach  we  return  to  a  direct  calculation  of 
N(E).  Because  we  also  calculate  D(E),  we  can  compare  the  two  approaches  and  judge  the  worth 
of  the  damage  energy  approximation  for  estimating  displacement  damage. 

In  the  next  two  sections  we  lay  out  the  formalism  that  leads  to  integral  equations  for 
Nv(E),  H(E),  B(E),  1(E),  which  we  have  discussed,  and  N|t(E),  the  number  of  replacements 
produced  by  an  ion  of  energy  E.  We  compare  our  results  first  with  the  binary-collision 
simulation  code  MARLOWE,^  and  then  with  the  standard  analytical  transport  theory  of  LNST. 

All  of  the  numerical  results  in  this  paper  and  much  of  the  disctission  assume  a  monoatomic 
crystal;  the  extention  to  complex  substrates  is  straightforward,  but  does  not  affect  the  points  made 
here. 

An  preliminary  version  of  this  work  was  reported  some  time  ago.’®  At  about  the  same 
time.  Coulter  and  Parkin’^  independently  investigated  some  of  these  same  questions. 

CROSS  SECTIONS  AND  STOPPING  POWERS;  GENERAL  FORMS 

Let  the  differential  cross  section  ♦(E,T)dT  represent  all  of  the  scattering  processes  that  an 
ion  can  undergo,  wherein  the  ion  retains  energy  E-T  and  transfers  energy  T  to  some  other 
particle  in  the  solid.  The  total  scattering  cross  section  is  given  by 
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i: 


♦(E)-  dT^(E.T) 


We  suppose  that  ♦(E.T)  is  composed  of  four  additive  parts. 


♦(E,T)  -  ♦e(E,T)  +  ♦„(E,T) 


♦e(E,T)  +  ♦t(E,T)  +  ♦b(E,T)  +  ♦v(E,T)  , 


(4 


(5 


representing  the  following  types  of  collisions:  electronic,  thermal,  replacement,  and  vacancy, 
respectively.  The  last  three  are  all  portions  of  the  elastic  scattering  cross  section  ♦m(£,T),  as 
pictured  in  fig.  1. 

The  electronic  cross  section  is  ♦E(E,T)dT,  where  T  represents  the  sum  of  the  energies 
transferred  from  the  moving  ion  to  the  electrons  of  the  lattice  atom.  AH  of  this  energy  goes  into 
inelastic  losses;  no  displacements  are  produced.  In  the  rest  of  this  paper  we  will  have  need  only 
for  the  electronic  stopping  power,  given  by 


Se(E) 


dT  T  ♦e(E,T)  , 


and  not  for  the  cross  section  itself. 

For  the  thermal  losses  we  write 


(6 


♦x(E,T)  dT  -  ♦h(E,T)  U(Eo-T)  dT  , 


(7 


where  U(  )  is  the  unit  step  function  and  ♦m  is  the  macroscopic,  elastic  scattering,  differential  cross 
section.  We  define  a  quantity  that  we  call  the  thermal  stopping  power. 


St(E) 


■E 

dT  T  ♦n(E,T)  U(Eo-T)  , 

0 


(8 


which  is  analogous  to  the  electronic  stopping  power. 

The  replacement  and  vacancy  cross  sections  are 


♦R(E,T)dT-  ♦n(E,T)  U(Eo-{E-T))  U(T-Eo)  dT 


(9 


and 

♦v(E,T)dT-  ♦m(E,T)  U(T-Eo)  U(E-T-Eo)  dT  .  (10 
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The  total  cross  sections  for  the  replacement,  vacancy  and  thermal  processes  are  given  by 


■i 


*,(E)-  dT*,(E,T),  i-s.v.r; 


and  we  set 


(11 


*rv(E) -♦»(£)  + *v(E)  •  (12 

In  the  same  spirit  that  we  defined  the  electronic  stopping  power  and  the  thermal  stopping 
power,  we  define  a  binding  energy  stopping  power,  given  by 


Sb(E)  -  Eg  ♦v(E)  . 


(13 


and  a  replacement  stopping  power,  which  represents  the  kinetic  energy  retained  by  the 
replacements  as  they  are  bound  onto  lattice  sites,  given  by 


S«(E) 


'E 

dT  (E-T)  «*(E,T)  . 

Jo 


(14 


DERIVATION  OF  THE  INTEGRAL  EQUATIONS 

The  following  discussion  parallels  the  approach  of  LNST.^^  Suppose  that  we  are 
interested  in  some  physical  process,  such  as  the  average  number  of  vacancies  Nv(E)  that  an  ion  of 
energy  E  will  produce.  Suppose  the  ion  undergoes  a  collision  that  results  in  its  retaining  an 
energy  E-T  and  transferring  energy  T  to  some  other  particle  in  the  solid.  These  two  particles  will 
traverse  the  solid  producing  more  vacancies;  we  measure  the  extent  of  these  effects  by  Nv(E-T) 
and  Ny(T-E||),  where  EK>£g  when  a  vacancy  is  produced  and  E|t>0  otherwise.  In  addition,  a 
vacancy  may  be  induced  at  the  collision  site  itself.  We  measure  this  collision  dependent  effect 
with  the  function  Nv(E,T),  which  will  be  either  zero  or  one. 

Because  some  such  collision  is  part  of  the  natural  chain  of  events  that  leads  our  original 
particle  to  produce  a  number  Nv(E)  of  vacancies,  we  know  that,  on  the  average,  the  effect  of 
Nv(E-T),  Nv(T-Eb)  and  Nv(E,T)  must  be  the  same  as  that  of  the  incoming  ion.  By  averaging  over 
the  collisional  possibilities,  we  can  express  this  equality  as 


«(E)  Nv(E)  - 


•E 

dT  *(E,T)  [  Nv(E-T)  +  Nv(T-Eb)  +  Nv(E,T)  ]  . 

Jo 


(15 


One  can  find  more  detailed  derivations  of  eq.  (!S)  in  LNST  and  elsewhere  (Sigmund  et  al,^^ 
Winterbon,  Sigmund  and  Sanders^^  (WSS)),  but  in  substance  they  reduce  to  this  balancing 
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argument 

If  we  split  the  cross  section  into  its  four  components  and  replace  the  total  cross  section 
4(E)  by  its  integral  form,  we  obtain 


E-Eo 

E-Eo 
Eo 

Eo 


dT  #*(E,T)  [  Nv(E)  -  Nv(T)  -  Nv.,(E,T)  ] 


dT  4v(E.T)  [  Ny(E)  -  Nv(E-T)  -  Nv(T-Eb)  -  Nv.v(E,T)  ] 


•1 

dT  4t(E,T)  (  Nv(E)  -  Nv(E-T)  -  Nv.t(E,T)  J 

.  0 


i: 


dT  4e(E,T)  (  Nv(E)  -  Nv(E-T)  -  Ny^^CE.T)  ]  , 


(16 


where  we  have  let  Nvj(E,T)  represent  the  number  of  vancancies  produced  at  the  collision  site 
when  the  j  (jav,R,T,E)  cross  section  is  operative.  Only  a  vacancy  producing  collision  results  in  a 
non-zero  value  of  Nyj(E,T),  j«v,R,T.E,  so  that 


Ny,y(E,T)  -  1 
Nyj(E,T)  -  0  ,  j««.T,E. 


(17 


In  eq.  (16)  LNST  would  have  replaced  Ny(T-Eg)  by  Ny(T),  because  Eg  is  small  compared 
to  the  other  energies  and  because  they  were  not  interested  in  the  effects  of  a  binding  energy.  We 
will  keep  the  Eg  parameter,  although  we  often  set  its  value  to  zero.  We  now  make  use  of  an 
approximation  of  LNST:  that  the  energy  transfer  in  electronic  collsions  is  small  compared  to  the 
energy  of  the  scattered  ion.  For  Ny(E-T)  they  would  write 

Ny(E-T)  »  Ny(E)  -  Nv'(E)  T  (18 

in  the  integral  with  the  electronic  cross  section,  so  that  the  first  two  integrals  in  the  last  line  of  eq. 
(16)  are  replaced  by  Sg(E)  Ny'(E).  We  also  impose  approximation  (18)  on  the  thermal  collisions,  so 
that  the  integral  of  the  first  two  terms  in  the  third  line  of  eq.  (16)  is  replaced  by  St(E)  Ny'(E). 
Because  the  energy  transferred  in  a  thermal  scattisring  is,  by  definition,  less  than  Eg,  the  relation 
(18)  will  be  a  good  approximation  in  thermal  collisions  except,  possibly,  for  energies  E<Eo.  We 
show,  in  Appendix  B,  that  this  approximation  1:  warranted.  For  convenience,  we  call  these  the 
electronic  stopping  power  and  thermal  stopping  power  approximations. 

When  we  make  use  of  the  various  definitions  of  cross  sections  and  stopping  powers  that 
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appear  in  the  previous  section  and  adopt  these  two  approximations,  we  can  write  our  integral 
equation  in  the  form 

(  Se(E)  +  St(E)  1  Nv'(E)  +  *,»(E)Nv(E)  +  J^E)  -  Gv(E)  (19 

where  Gv(E)  is  the  integral  expression 

Gv(E)  -  dT  {  ♦«(E,T)  Nv(T)  +  ♦v(E,T)  (  Nv(E-T)  +  NylT-Eg)  ]  )  .  (20 

and  Jv(E)  is  the  inhomogeneous  term 

Jv(E)-  -♦v(E).  (21 

In  exactly  the  same  manner  (see  Appendix  A)  we  can  derive  integral  equations  for  the 
quantities  Nk,  B,  H,  and  I.  We  solve  these  integral  equations  numerically  using  a  cubic  spline 
method 

COMPARISON  WITH  BINARY-COLLISION  SIMULATIONS 

The  best,  readily  available,  method  of  calculating  displacement  damage  is  the  binary- 
collision  simulation  code  MARLOWE,^  which  simulates  radiation  effects  by  following  individual 
ions  through  a  psuedo  crystal.  Among  the  many  differences  between  the  simulation  code 
MARLOWE  and  the  present,  analytical  transport  theory,  approach  are  the  following: 

a)  MARLOWE  modeb  a  crystalline  material,  instead  of  treating  the  substrate  as 
amorphous. 

b)  MARLOWE  calculates  the  exact  scattering  angles  and  energy  transfers  for  each 
collbion  using  the  Moliere  potential.^  Thb  potential  is  a  version  ofthe  Thomas- 
Fermi  potential  with  an  exponential  tail  and  is  thought  to  be  more  realistic.^^  The 
present  theory  uses  an  approximate  cross  section^  based  on  the  Moliere  potential, 
which  b  obtained  by  using  the  same  approximations^  that  Lindhard,  Nieben  and 
Scharff^^  (LNS)  use  to  obtain  their  cross  section  from  the  Thomas- Fermi  potential. 
The  detailed  forms  of  the  Moliere  potential  and  cross  section  can  be  found  in 
Appendix  C. 

c)  With  MARLOWE  one  has  a  choice  of  inelastic  loss  mechanisms.  Several 
differential  cross  sections  for  electronic  loss  are  built  into  the  code,  as  is  a 
nonlocal,  continuous  slowing  down,  model  of  electronic  less.  The  present 
analytical  method  incorporates  only  a  continuous  slowing  lown  model  that  is 
similar,  but  not  identical  to,  that  used  in  MARLOWE. 

d)  MARLOWE  imitates  the  effects  of  many  body  collisions  by  allowing,  in  a  fashion, 
several  binary  collbions  to  occur  simultaneously. 
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e)  MARLOWE  imposes  a  cutoff  in  the  maximum  impact  parameter  for  which  a 
collision,  elastic  or  inelastic,  may  occur. 

f)  MARLOWE  remembers  previous  damage,  so  that,  if  the  cascade  doubles  back  on 
inseif,  the  same  lattice  atom  cannot  be  removed  from  its  lattice  site  twice.  The 
analytical  code  is  essentially  an  ensemble  average,  each  (figurative)  element  of 
which  sees  only  virgin  substrate. 

As  the  first  stage  of  comparison,  we  want  to  eliminate  as  many  of  these  differences  as 
possible.  In  particular,  if  we  could  turn  all  of  them  off,  and  if  we  then  found  that  the  two 
calculations  gave  similar  results,  we  could  examine  the  effect  of  each  dissimilarity  by  selectively 
turning  it  back  on. 

We  examine  the  case  of  an  fee  copper  crystal.  For  the  screening  length  that  appears  in  the 
Molidre  potential  we  use  the  value  a>7.38  pm.^  For  the  energy  parameters  discussed  in  the 
introduction  we  use  Eg^o  and  £^>£0^25  eV.  In  order  to  effect  the  comparison,  we  make  the 
following  modifications  and  adjustments  to  MARLOWE: 

a)  We  run  MARLOWE  in  its  amorphous  mode,  in  which  the  crystal  is  rotated  through 
some  random  angle  after  each  scattering.  This  has  the  effect  of  imitating,  if  not 
duplicating,  the  amorphous  character  of  the  analytical  code. 

b)  We  insert  into  MARLOWE  a  new  scattering  angle  subroutine  that  simulates  the 
effect  of  the  approximate  cross  section  that  the  analytical  code  uses.  In 
Appendix  C  we  indicate  how  to  start  with  the  approximate  Moliere  cross  section 
and  obtain  the  scattering  angle  versus  impact  parameter  relation  that  we  use  in 
MARLOWE. 

c)  We  turn  off  the  inelastic  stopping  in  both  codes.  Because  of  the  conceptual 
differences  in  the  two  models  of  inelastic  loss,  there  is  no  simple  way  to  make 
them  comparable.  In  any  case,  we  are  not  interested  at  the  moment  in  differences 
in  displacement  damage  caused  by  different  electronic  loss  models. 

d)  We  make  use  of  a  MARLOWE  provision  that  allows  the  simultaneous  collision 
feature  to  be  turned  off. 

To  a  greater  or  lesser  extent  these  modifications  remove  the  major  dissimilarities  between  the  two 
approaches. 

For  the  moment  we  ignore  differences  (e)  and  (f),  as  well  as  any  others  that  we  have  not 
identified.  In  order  to  compare  the  two  approaches  we  plot  the  displacement  efficiency  as  a 
function  of  energy.  The  displacement  efficiency  K^^E)  is  defined  in  terms  of  the  total  number  of 
displacements  by 

N(E)  -  Kv(E)  EKE)  /  2Eo  .  ( '.2 

The  quantity  Kv(E)  measures  the  efficiency  with  which  vacancies  are  produced  relative  to  the 
hard  core  model,  for  which  Ky(E)-l.  Because  we  have  turned  off  the  inelastic  losses,  all  of  the 
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DISPLACEMENT  EFFICIENCY  [Kv|E)J 


energy  is  damage  energy,  so  that  D(E)-E. 

In  fig.  2  we  show  that  displacement  efficiency  for  the  present  theory,  for  MARLOWE 
with  the  changes  listed  above  and  for  a  straight  MARLOWE  calculation  modified  only  by  the 
absence  of  inelastic  losses.  We  see  that  the  agreement  between  the  present  theory  and  the 
modified  MARLOWE  results  is  quite  good.  Both  approaches  indicate  that  the  displacement 
efficiency  becomes  constant  at  0.5- 1.0  keV.  For  the  standard  MARLOWE  calculation  there  is 
more  structure  in  the  displacement  efficiency  curve  due  to  the  effects  of  the  lattice.  The  more 
realistic  MARLOWE  calculation  indicates  a  somewhat  (9  percent)  higher  disolacement  efficiency 
above  1  keV. 

To  study  the  exact  origin  of  this  difference,  we  start  with  the  standard  MARLOWE 


ENERGY  IE  (keV)l 


Fig.  2  —  A  comparison  of  Jisplacement  efficiency  as  estimated  by  the  present  theory, 
MARLOWE,  and  MARLOWE  as  modified  to  contain  some  of  the  same  approximations  as  the 
present  theory.  (The  case  shown  is  for  fee  copper  with  Eds25  eV  and  no  inelastic  losses.) 
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calculation,  at  £»!  kev,  and  selectively  turn  on  various  combinations  of  approximations.  Table  1 
shows  the  results  of  this  procedure.  We  see  that  the  effects  of  the  different  approximations  are 
not  additive.  The  only  consistent  effect  seems  to  be  that  an  amorphous  substrate  tends  to  reduce 
the  number  of  vacancies  by  S-10  percent. 


Table  1 

—  Effect  of  Various  Approximations  on  the 
Number  of  Vacancies 

Approximations  Added  to  the  Standard 

Percentage  Change 

MARLOWE  Calculation 

in  N(E)  at  1  keV 

Amorphous 

Isolated 

WSS  Cross 

Substrate 

Collision 

Section 

- 

X 

- 

.  +1 

- 

- 

X 

.  +3 

- 

X 

X 

.  0 

X 

- 

- 

.  -4 

X 

X 

- 

.  -5 

X 

- 

X 

.  -11 

X 

X 

X 

.  -9 

We  can  compare  the  efficiency  with  which  replacements  are  produced  in  the  same  manner 
in  which  we  examined  vacancies.  We  define  a  replacement  efficiency 

Nr{E)  -  Kr(E)  EKE)  /  2Eo  .  (23 

Using  the  same  assumptions,  we  again  have  0(E)»E.  We  see  in  fig.  3  that  the  present  results  are 
somewhat  lower  but  of  the  same  pattern  as  the  amorphous,  WSS  cross  section,  isolated  collision, 
version  of  MARLOWE.  But  the  effect  of  neglecting  the  crystal  structure  is  large  in  this  case,  as 
one  might  expect;  the  present  results  are  only  30  percent  of  the  straight  MARLOWE  results. 


11 


REPLACEMENT  EFFICIENCY  IK, (Ell 


ENERGY  (E  (keV)] 


Fig.  3  —  A  comparison  of  replacement  efficiency  shown  for  the  same  conditions  as  the  last 

figure. 


DAMAGE  FRACTION  (D(E|/E] 


COMPARISON  OF  DAMAGE  ENERGY  ESTIMATES 

Comparisons  of  the  present  approach  and  MARLOWE  are  difficult  when  electronic  losses 
are  included.  We  use  the  LNST  continuous  slowing  down  model  (see  eq.  (CIS));  the  Firsov 
stopping  power  coefficient  has  the  value  ev^^^  A'^  for  copper.  The  energy  loss  process 

is  applied  continuously  in  the  differential  equation  so  that,  figuratively,  each  moving  ion 
continuously  undergoes  inelastic  losses.  MARLOWE  can  also  use  this  non-local  stopping  law,  but 
it  is  applied  to  the  line  segments  that  are  generated  by  the  elastic  collision  asymptotes.  At  low 
energies,  especially,  there  are  a  number  of  differences  in  the  two  methods. 

MARLOWE  also  can  use  a  local  electronic  loss  mechanism,  with  a  differential  cross 
section  of  exponential  form,  such  that  the  stopping  power  resulting  from  that  cross  section  is 
essentially  the  same  as  that  of  the  continuous  slowing  down  model.  In  fig.  4  we  show  the  results 
of  MARLOWE  runs  with  both  stopping  mechanisms  and  the  results  of  the  present  approach.  At 
low  energies  there  is  notable  disagreement  between  them.  At  higher  energies  the  local  stopping 


0.03  0.1  0.3  1.0  3.0  10 


ENERGY  [E  (keV)l 

Fig.  4  —  A  comparsiou  of  the  damage  fraction  (D(E)/E)  for  the  same  conditions  as  fig.  2, 

except  that  inelastic  losses  are  included. 
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DAMAGE  FRACTION  (D(E)/E) 


method  in  MARLOWE  agrees  more  closely  with  the  present  approach  than  with  the  other 
MARLOWE  model.  We  conclude  that  differences  in  the  models  preclude  good  agreement  in 
estimates  of  damage  energy. 

EFFECT  OF  POTENTIAL 

The  conventional  LNST  energy  partitioning  theory  is  based  on  the  Thomas- Fermi 
interaction.  Figure  5  serves  to  demonstrate  the  effect  on  energy  partitioning  of  changing  the 
interaction.  All  of  the  curves  are  based  on  the  formalism  of  this  paper.  The  solid  curve 
represents  the  results  of  using  the  WSS’*  fit  to  the  LNS^*  version  of  the  Thomas-Fermi  cross 
section,  with  the  Lindhard  screening  length  a- 10.78  pm.  (The  details  of  these  interactions  and  the 
expressions  for  the  screening  lengths  are  given  in  Appendix  C.)  The  two  dashed  curves  represent 


ENERGY  [E  (keV)] 

Fig.  5  —  A  comparison,  using  the  present  theory,  of  the  damage  fraction  ic  copper  using  three 

different  cross  sections. 


the  results  of  using  the  cross  section  based  on  the  Moliire  potential,  both  for  the  Lindhard  value 
of  the  screening  length  and  the  Torrens  and  Robinson^  value  a-7.38  pm,  which  was  obtained  by 
adjusting  the  screening  length  so  that  the  low  energy  portion  of  the  Moliere  potential  agreed  with 
a  Cu-Cu  interaction  that  reproduced  some  of  the  bulk  properties  of  copper.  We  see  that  there  are 
large  shifts  away  from  the  LNST  (Thomas-Fermi  interaction  with  LNST  screening  length)  in  the 
damage  energy,  produced  by  switching  from  the  Thomas-Fermi  to  the  more  realistic  Moliere 
potential,  and  then  produced  by  switching  to  the  smaller,  more  realistic,  screening  length. 

We  can  make  another  comparison  of  the  effects  of  changing  the  potential  be  examining 
the  number  of  vacancies  produced  in  each  case.  In  fig.  6  we  plot  the  number  of  vacancies 
produced  per  unit  energy  (N(E)/E)  as  a  function  of  energy.  The  two  potentials  that  use  the  10.78 
pm  screening  length  are  in  rough  agreement  at  all  energies,  presumably  because  the  Thomas- 
Fermi  and  Moliere  cross  sections  coincide  for  higher  relative  energy  transfers.  Upon  comparing 


Fig.  6  —  A  comparison  of  the  number  of  vacancies  produced  for  the  same  three  cross  sections 

as  the  last  figure. 


fig.  S  and  fig.  6,  we  fail  to  see  any  clear  relation  between  damage  energy  and  number  of 
vacancies,  relative  to  changes  in  the  potential. 

We  can  examine  this  disparity  by  looking  at  other  properties  of  the  cross  sections.  In 
fig.  7  we  show  the  ratio  of  the  total  elastic  stopping  power  to  the  electronic  stopping  power  for 
each  of  the  three  potentials.  The  total  elastic  stopping  is  given  by 


Sa(E) 


dT  T  *m(E,T)  . 


(24 


Clearly,  the  reason  that  the  Thomas-Fermi  interaction  leads  to  greater  damage  energy  than  does 
the  Moliere  is  that  the  elastic  stopping  dominates  the  electronic  to  a  much  greater  extent  in  the 
Thomas-Fermi  case.  We  note  that  the  WSS  cross  section  and  the  Moliere  (each  with  the  10.78  pm 
screening  length)  lead  to  the  same  stopping  at  high  energies,  as  we  would  expect.  But  the 
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Fig.  7  —  A  comparison  of  the  ratio  of  nuclear  and  inelastic  stopping  powers  for  the  same  three 

cross  sections  as  fig.  S. 
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PERCENTAGE  CHANGE 


corresponding  damage  energy  curves  do  not  merge.  The  relative  amounts  of  elastic  and  inelastic 
loss  depends  sensitively  on  the  behavior  of  the  many  cascade  atoms  moving  with  lower  energies, 
for  which  the  disparity  in  relative  stopping  power  is  still  felt. 


BINDING  AND  CUTOFF  EFFECTS 

In  fig.  8  we  show  the  effect  of  introducing  a  binding  energy  of  10  eV  on  the  various 
measures  of  damage.  (This  is  an  unnaturally  large  value,  but  it  serves  to  indicate  tendencies.)  We 
assume  that  each  displacement,  but  not  replacement,  costs  the  displaced  atom  10  eV.  The  solid 
curve  in  fig.  8  shows  the  percentage  reduction  in  the  number  of  Frenkel  pairs  as  a  function  of 
energy.  Abo  shown  are  the  percentage  changes  in  the  thermal  loss  H(E)  and  the  damage  energy 
D(E)>H(E)-»-B(E).  We  see  that  the  imposition  of  a  binding  energy  produces  opposite  effects  in  the 
damage  energy  and  the  number  of  dbplacements,  indicating  that  damage  energy  is  not  a  good 
measure  of  dbplacement  production. 

The  imposition  of  a  finite  cutoff  energy  makes  no  change  in  the  number  of  Frenkel  pairs 
produced,  but  the  asymptotic  value  of  the  damage  energy  fraction  b  reduced  by  18  percent  and 
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Fig.  8  —  The  percentage  change  in  the  number  of  vacancies,  damage  energy  and  thermal  losses 

when  a  10  eV  binding  energy  is  imposed. 
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25  percent  when  the  cutoff  energy  is  25  eV  and  50  eV,  respectively.  This  is  another  example  that 
the  estimates  of  the  number  of  vacancies  and  of  the  damage  energy  are  sensitive  to  different 
effects. 

COMPARISON  TO  LNST 

Our  formalism  coincides  with  that  of  LNST  when  E^,  Eg  and  Eg  are  zero,  in  which  case 
our  integral  equation  for  D(R)  is  identical  to  the  LNST  equation  for  damage  energy.  In  the 
absence  of  a  finite  displacement  energy,  we  would  naively  expect  that  there  would  be  no  damage 
energy,  because  all  of  an  ion’s  energy  would  eventually  be  drawn  off  by  collisions  with  the 
electrons.  LNST  find,  however,  that  if  the  inelastic  stopping  goes  to  zero  as  E^  and  if  the 
interatomic  potential  falls  off  as  r'*,  then  there  is  a  non-zero  solution  of  the  damage  energy 
equation  provided  that  s<2/(l-p).  (In  fact,  the  value  of  s  should  not  be  very  close  to  the  limit.) 
For  the  standard  velocity-proportional  inelastic  stopping,  this  requires  s<4.  LNST  use  the 
Lindhard,  Nielsen  and  Scharff  cross  section,  which  is  a  form  of  the  Thomas- Fermi  cross  section 
that  has  the  low  energy  behavior  equivalent  to  s>3.  In  order  to  effect  a  comparison  with  LNST, 
we  use  the  WSS  fit  to  the  LNS  cross  section  (see  Appendix  C,  eq.  (Cl)  and  eq.  (C8))  and  the  LNST 
electronic  stopping  (eqs.  (C15-C17)). 

In  fig,  9  we  compare  our  results  for  the  fraction  of  energy  that  is  damage  energy.  The 
LNST  curve  is  the  equal  mass  result  of  Winterbon,^  plotted  as  a  function  of  the  reduced  energy 
value  «-E/El  (Eq.  (C4)).  The  present  theory  is  represented  by  three  curves,  for  three  different 
values  of  «o-Eo/El.  (The  limited  extent  of  the  €o»5xlO*^  curve  is  due  to  a  numerical  quirk  in 
the  code  that  was  not  worth  removing  for  the  sake  of  this  one  demonstration.)  For  definiteness, 
in  the  case  of  copper  the  specific  values  of  the  displacement  energy  corresponding  to  the  the  three 
curves  are  0.14  eV,  14  eV  and  1400  eV.  Even  in  the  case  of  the  tiny  0.14  eV  displacement 
energy,  there  is  a  marked  difference  between  our  results  and  the  Eg-O  results  of  the  LNST 
method. 

The  LNST  results  are  always  noticeably  lower  than  the  present  results.  The  reason  for  this 
is  that  in  our  approach,  for  energies  below  Eg,  certain  modes  of  energy  loss  are  prohibited.  Most 
ions  with  energies  below  Eg  are  tied  to  lattice  sites,  either  by  way  of  having  been  a  lattice  atom 
that  received  a  subthreshold  energy  transfer  or  because  they  arrived  at  a  lattice  site  as  a 
replacement.  In  both  cases,  in  our  model,  they  are  incapable  of  losing  energy  to  electronic 
processes.  In  the  LNST  model  ions  continue  to  lose  energy  electronically  until  they  have  zero 
energy.  This  effect  is  insensitive  to  changes  in  the  displacement  energy.  When  we  halve  the 
displacement  energy,  we  roughly  double  the  number  of  displaced  atoms,  each  of  whose  energy 
eventtially  drops  below  the  new  Eg  value  and  is  prohibited  from  losing  energy  electronically.  It  is 
this  prohibit  k1  loss  mechanism  in  the  present  approach  that  accounts  for  the  qualitative  difference 
between  the  curves  in  fig.  9. 

In  closing  this  section,  we  note  that  there  are  other  ways  to  estimate  a  damage  energy 
without  imposing  restrictions  on  the  low  energy  behavior  of  the  interatomic  potential.  Miller  and 
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Boring,^  for  instance,  require  that  the  inelastic  stopping  vanish  at  low  energies,  of  the  the  order 
of  the  ionization  energies  of  the  electrons.  As  in  our  model,  this  condition  prohibits  inelastic 
losses  below  certain  energies  and  guarantees  a  finite  damage  energy. 


CONCLUSIONS 

The  present  theory  yields  good  estimates  of  displacement  damage  when  compared  to  the 
binary-collision  simulation  code  MARLOWE.  Agreement  is  quite  good  when  MARLOWE  is  run 
in  its  amorphous  mode  using  the  modified  cross  section  that  is  used  by  the  present  approach.  But 
even  when  these  allowances  are  not  made,  the  present  theory  differs  from  MARLOWE  by  less 
than  10  percent  in  the  case  (fee  copper)  we  examined.  Changes  in  interatomic  potential  have  a 
much  different  effect  on  the  damage  energy,  which  is  widely  used  to  estimate  displacement 
damage,  than  they  do  on  a  direct  calculation  of  the  number  of  displacements.  Similarly,  a  finite 
binding  energy  creates  differing  responses  in  the  two  approaches. 


REDUCED  ENERGY  U  »  E/EJ 

Fig.  9  —  A  comparison  of  the  damage  fraction  as  a  function  of  energy  for  a  variety  of 

displacement  energies. 


Finally,  we  have  noted  that  any  evaluation  of  the  damage  energy  requires  some  restriction 
on  the  relative  strength  of  the  inelastic  and  elastic  cross  sections  at  low  energies.  LNST  assume  a 
very  soft  low-energy  elastic  cross  section,  whereas  we  assume  that  atoms  bound  to  the  lattice 
suffer  no  inelastic  loss.  We  again  find  that  these  differences  in  the  models  create  large 
differences  in  the  estimates  of  damage  energy. 

We  conclude  that  the  direct  calculation  of  the  displacement  production  is  considerably 
more  reliable  than  estimating  it  from  the  damage  energy. 
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APPENDIX  A;  DERIVATION  OF  THE  OTHER  INTEGRAL  EQUATIONS 

In  the  same  manner  that  the  integral  equation  for  Nv(E)  was  derived  in  Section  3,  we  can 
develop  the  equivalent  integral  equations  for  the  other  quantities  of  interest.  We  can  let  Fj(E) 
(i-R,B,T,E)  stand  for  the  quantities  N*(E),  B(E),  H(E),  and  1(E),  repsectively.  Then  F,-^j(E,T) 
represents  the  residual  quantity  F;(E,T)  that  appears  when  a  collision  occurs  under  the  influence 
of  the  j  cross  section,  j«v,R,T,E.  The  only  non-zero  F,-  j(E,T)  are 

Fv,v  -  1 
Fr.r  -  1 
Fr.v  -  Eb 
Fr.r  -  T 
Ft.r  -  E-T 

Fe.e-T.  (A1 

The  form  (19)  can  be  used  to  represent  the  integral  equation  for  F{(e), 

(  Se(E)  +  St(E)  ]  Fi'(E)  +  ♦«,(£)  Fi(E)  +  J,(E)  -  Gi(E)  (A2 

where  Gf(E)  is  the  integral  expression 

G,(E)  -  I  ^  dT  {  ♦r(E,T)  F,(T)  +  #v(E,T)  [  Fj(E-T)  +  F^T-Eg)  ]  }  .  (A3 

Jo 

and  the  appropiate  inhomogeneous  term  for  these  are 

Ji(E)-  -*r(E),  -Sb(E),  -St(E)-Sr(E),  -Se(E),  (A4 

for  i-R,B,T,E.  We  note  that  if  the  equations  for  B(E),  H(E)  and  1(E)  are  added,  we  obtain  an 
integral  equation  for  F(E)-B(E)+H(E)+I(E)  that  is  satisfied  by  F(E)»E. 
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APPENDIX  B:  TEST  OF  THE  THERMAL  STOPPING  APPROXIMATION 

In  order  to  assess  the  accuracy  of  the  thermal  stopping  approximation  we  consider  its 
effect  on  the  calculation  of  1(E).  Replacing  the  integral 

fEo 

dT  ♦t(E,T)  I  1(E)  -  I(E-T)  ] 

Jo 

by  St(E)  r(E)  will  clearly  be  least  accurate  for  E<Eo.  For  energies  of  this  order  we  can  easily 
obtain  the  first  few  terms  of  a  power  series  solution^  of  the  integral  equation  for  1(E).  Using  the 
WSS  (Winterbon,  Sigmund  and  Sanders)^*  form  of  the  LNS^^  version  of  the  Thomas-Fermi  cross 
section  (see  Appendix  C,  Eq.  (Cl)  and  Eq.  (C8)),  and  using  the  LNST^^  electronic  stopping  (Eqs. 
(CIS- 16)),  both  solutions  can  be  written  (in  the  equal  mass  case)  in  the  form 

1(E)  /  E  -  ai  s  +  a2  s^  +  as  s^, 

s  -  (2k/A)  (E/Ej’''*  (B1 

where  Ey.  is  the  Lindhard  energy  unit  (eq.  (C4)),  k  is  the  electronic  stopping  coefficient 
(Eq.  (CIS)),  and  A*  1.309  appears  in  the  WSS  cross  section.  We  find  that  the  coefficients  (a^  a?, 
as)  are  (0.S973,  -0.3794,  0.2S43)  for  the  present  theory  without  the  thermal  stopping 
approximation  and  are  (0.S7 14,-0.3333,  0.197S)  with  the  approximation.  The  form  (Bl)  should 
only  be  used  for  energies  less  than  the  displacement  energy,  due  to  the  discontinuity  in  the  second 
derivative  of  1(E)  at  E«Eo. 

When  we  use  typical  values  of  Eg  and  examine  elements  from  beryllium  to  uranium,  we 
find  a  systematic  error  in  1(E)  of  about  4  percent  at  E-Eg  due  to  the  thermal  stopping 
approximation.  The  error  for  energies  less  than  Eg  is  only  slightly  larger.  It  should  be  negligible 
at  higher  energies.  This  error  is  certainly  less  than  other  inaccuracies  associated  with  this  type  of 
displacement  damage  calculation,  so  we  are  justified  in  using  the  thermal  stopping  approximation. 
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APPENDIX  C;  SPECmC  CROSS  SECTIONS  AND  POTENTIALS 

The  most  commonly  used  cross  section  for  atomic  scattering  problems  is  the  LNS^^  version 
of  the  Thomas-Fermi  cross  section.  By  using  what  we  will  call  the  LNS  method,^  they  obtain  a 
simple,  universal,  cross  section  based  on  the  Thomas-Fermi  potential.  They  write  the  cross 


section  as 

♦h(E,T)  dT  -  n  »  3-2  f^(3)  ^j3  ^  (C, 

s2  -  E  T  /  (  7  El^  )  ,  (C2 

T  4  A,  A2/  (Ai+A2)2  ,  (C3 

■  Zi  Z2  (  Ai  +  A2  )  /  (  n  A2  )  (C4 

a  -  [2(3x/32)2/^]  Z*’'^  «  0.8853  ao  Z*'''^  ,  (C5 

Z  -  Zl  -  (  Zi^/^  +  22^'^  .  (C6 

where  n  is  the  number  density  of  the  atoms  in  the  substrate.  The  LNS  screening  length  depends 
on  their  choice  of  Z>Z|..  Another  common  choice  is  that  of  Firsov,^ 

Z  -  2f  -  (  2/^2  ^  (C7 

LNS  supply  the  kernel  of  their  cross  section,  f^Is),  in  tabular  form.  More  convenient  is  the  fit  to 
fj.  by  WSS,  in  the  form’® 

fg(s)  -  A  s’-'®  (  1  +  (2  A  ,  A  »  1.309  .  (C8 

The  Moliere  potential  is^® 

V(r)  -  I  Z,  Z2  e^  /  r  1  x«(r/a)  .  (C9 

Xk(x)  -  0.35  exp(-0.3x)  +  0.55  exp(-1.2x)  +  0.10  exp(-6x)  (CIO 


One  can  follow  the  LNS  method  and  obtain  the  cross  section  equivalent  (within  the  assumptions 
of  the  LNS  method)  to  this  potential.  This  cross  section  has  been  fit  to  a  convenient  form^^ 

f  (0.007  +  s/2)  /  (0.0387  +  0.826  s  ♦  s^  )  ,  s>s*-.06 
fH(s)-j  ,  3  ^  (Cll 

[  -s  20.45  ln(s)  -  71  s  +  422.097  s^  -  1429.70  s^  ,  s<s* 
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In  order  to  use  the  LNS  approximate  form  of  the  Moli&re  cross  section  in  the  code 
MARLOWE,  we  make  use  of  the  relation  between  the  impact  parameter  and  the  energy  transfer 
given  by 


(p/a)*  - 


'oo 

ds  s'*  f(s) 

.  s* 


(C12 


where 


s*  -  (E  T*  /  7  El*)’'*  .  (C13 

Given  the  value  of  s*  in  eq.  (Cl 2)  corresponding  to  p/a,  we  can  determine  the  energy  transferred 
T*  in  the  collision  from  eq.  (C13).  The  scattering  angle  is  then  given  by 


-  2  sin*’ (T*/7E)’'*  .  (Cl 4 

We  note  that  we  also  set  the  apsis  of  the  collision  equal  to  the  impact  parameter,  which  is  a  part 
of  the  LNS  approximation. 

The  LNST’*  expression  for  of  the  electronic  stopping  (eq.  (6))  is 

Sh(E)  -  kL  E’'*  ,  (Cl 5 

kL  -  .7323  n  Z/'*  Zj  Zl*’  Aj*’  A,'’'*  eV’'*  A*’  .  (C16 

where  n  is  in  units  of  gm/cm^  and  the  masses  in  amu.  The  equivalent  Firsov  expression  can  be 
obtained  from  Eqs.  (16-19)  of  Ref.  9;  it  also  has  the  form  (C15)  with 

kF-.2024  n  (  Z,  +  Z2)  Aj'’  Ai*’'*  eV’'*  A*’ .  (C17 
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